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Abstract 

We present the implementation of an algorithm for graph isomor- 
phism testing, based on ideas about number of walks (of sufficiently 
large length) between vertices, expanded for strongly regular graphs 
(SRG-s) by testing the local complements and values of determinants 
of their adjacency matrices. All known non-isomorphic SRG-s (with 
no more than 64 vertices) are distinguishable by this method. 
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1 Introduction, basic algorithm 

The graph isomorphism problem is one of difficult central problems in graph 
theory. It can be easily formulated as the question, whether two differently 
looking images represent the same graph. We will consider only undirected 
graphs without loops and multiple edges. Let us denote by V(G), E{G) the 
sets of vertices and edges of graph G. Two graphs G, H are isomorphic iff 
there exists a bijective mapping / : V(G) i— ► V(H) between sets of vertices 
of graphs G and H such that for all pairs of vertices u, v G V(G) the edge 
{it, v} belongs to E(G) if and only if {f(u), f(v)} is an edge of graph H. 

The graph invariant A is a mapping from the class of all graphs Q to 
the set 1Z such that for any isomorphic graphs G, H the images X(G), \(H) 
are the same, i.e. X(G) = \{H). Commonly used graph invariants are e.g. 
number of vertices, number of edges, degree-sequence (the non- descending 



sequence of degrees of the vertices), eigenvalues or characteristic polynomial 
of the adjacency matrix, etc. All above-mentioned invariants are incomplete, 
which means that we can find pairs of nonisomorphic graphs with the same 
invariants. 

There are also complete graph invariants, i.e. A(G) = X(H) holds if 
and only if graphs G, H are isomorphic (see e.g. [T] for greatest charac- 
teristic string invariant). Such invariants are not known to be computable 
in polynomial time, because the existence of polynomial algorithm for the 
class of graph isomorphism problems is still an open question. The existing 
polynomial algorithms are inexact in the sense that they cannot recognize 
some pairs of nonisomorphic graphs (because they use incomplete invariants 
for testing). Our proposed algorithm also belongs to this type of non-exact 
polynomial algorithms. 

Let us briefly summarize the main results from Czimmermann [Sj 
where the basic variant of isomorphism testing algorithm using walks between 
vertices counting was formulated. 

Let A = {dij) be the adjacency matrix of graph G. Then the number of 
all possible walks of length k from vertex v$ G V(G) to vertex Vj G V(G) is 
simply the element of matrix power A k . Let us assign to every vertex 
Vi G V(G) the set 



i.e. Sij is the sequence formed from number of walks of length 1,2, ... ,m 
which begin in vertex Vi and end in Vj. Let us have a pair of graphs G,H 
with equal number n of vertices for which the sequences of sets (0) are 



If G, H are isomorphic, then for suitable permutation of vertices of graph 
G the two sequences are equal. The proposed algorithm is based on com- 
parison of permuted sequences Wg, V^h with m (maximal length of walks 
considered) sufficiently large. In [3] it was shown that m can be chosen as 
the number of common distinct eigenvalues of graphs G, H (there is no addi- 
tional information coming from walks of greater length). For this basic form 
of algorithm it is easy to find the causes when it fails; e.g. two nonisomorpic 
strongly regular graphs with the same set of parameters are indistinguishable. 
The graph G is called strongly regular (SRG) with parameters n, d, a, j3 if it 
has n vertices, every vertex has degree d, every pair of adjacent vertices has 
a common neighbours and every pair of distinct nonadjacent vertices has j3 
common neighbours. For example, the well-known Petersen graph (see Fig. 
HJ) is SRG with parameters (10,3,0,1). 






W G = {S 1 ,S 2 ,...S n }, W H = {T U T 2 ,..., 



T n }. 
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Figure 1: Petersen graph - strongly regular graph with parameters (10,3,0,1) 

Each SRG has only three distinct eigenvalues, which are completely de- 
termined by parameters n, d, a, (3. Thus, the complete information about 
walks is contained in sequences of sets of walks with length k < 3. Such 
short walks are not sufficient to discover non-isomorfism of SRG-s with the 
same set of parameters. For SRG-s, we need to investigate some additional 
properties. 

2 Improved algorithm 

for strongly regular graphs 

The key idea is to destroy the regularity of given SRG by representing it with 
the sequence of suitable modified non-regular graphs. We will explore the 
notion of local complement (see for this purpose. 

Definition 2.1 Let G be a graph. The local complement Gl{u) of G at the 
vertex u G V(G) is defined to be the graph with the same vertex set as G such 
that: 

1. if v and w are distinct neighbours ofu, then they are adjacent in Gl(u) 
if and only if they are not adjacent in G, 

2. if v and w are distinct vertices of G, and not both neighbours of u, then 
they are adjacent in Gl{u) if and only if they are adjacent in G. 

Now, instead of testing isomorphism of two (SRG) graphs G, H, we can 
test for existence of suitable permutation of vertices of graph H such that 
the sequences of local complements of G, H have the same graph invariants 
(e.g. the sequences of sets of walks, mentioned in previous section). 
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Table 1: Groups 
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nonisomorfic SRG-s 



In Table ^ there are parameters and numbers of nonisomorphic SRG- 
s with given parameter sets for up to 64 vertices. We obtained adjacency 
matrices for all mentioned graphs from WEB-sources j7j, [Sj- To make the 
later processing easy, the SQL database of SRG-s were created with the aid of 
simple but powerful database engine SQLite (see [9 ) and the corresponding 
interface module pysqlite ^0] for programming language Python. 

In the paper jH] the authors were able to distinguish all graphs from this 
table (although they have omitted two large groups of SRG-s with parameters 
(35,16,6,8) and (37,18,8,9)) with their algorithm based on eigenvalues of 
certain matrix inspired by the notion of coined quantum walks. We will 
show that our modified algorithm can distinguish all of the graphs in Table 
^ when combined with such easy to compute incomplete graph invariants as 
determinants of adjacency matrices of local complements. 



3 Algorithm implementation and testing 

We will focus on testing our algorithms only for set of SRG-s from Table ^ 
because strongly regular graphs cause main difficulties for graph isomorphism 
algorithms. 

For the rapid implementation of basic algorithm (see section 1) we choose 
the interpreted programming language Python [12] ■ It has very clear syntax 
and together with SciPy (Scientific tools for Python, ^3]) and IPython (en- 
hanced Interactive Python shell, see (T^j) it forms very powerful and conve- 
nient environment for rapid development of prototype applications. 

In the basic algorithm we need to compute the powers A k of 0-1 matrices 
A (for SRG-s, up to k = 3). Underlying C and FORTRAN libraries in SciPy 
use machine integer and floating point numbers with fixed precision. This 
will give wrong results when the elements of A k will be large (e.g. the 17-th 
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power of adjacency matrix of the unique SRG with parameters (36, 10,4, 2) 
shows on 32-bit computers some negative entries, which is completely wrong). 
The same problem exists with computing determinants (for the above matrix 
even on 64-bit computer we see the nonsensical value -351843720888319.81 
of determinant). 

Python supports arbitrary precision integers, but it has no efficient lin- 
ear algebra operations with integer matrices and vectors. So we decided to 
complement the software tools used with NTL j!5j . a C++ library for manip- 
ulating arbitrary length integers, and for vectors, matrices, and polynomials 
over such integers. We used the ctypes Python module (see [Tfi] ) to call the 
functions in NTL library in Python. 

We verified all "small groups" (with at most 227 members) from Table 
H using only basic algorithm applied to local complements. But the testing 
for three "big" groups of SRG-s was very time-consuming, so we considered 
to make preliminary classification of nonisomorphic graphs with the same 
parameters, using some cheap graph invariants. 

In Czimmermann j3] there is shown that if two graphs G, H have distinct 
determinants of their corresponding adjacency matrices, then the spectra of 
adjacency matrices are also distinct and that two graphs with distinct spectra 
can be distinguished by the basic variant of algorithm, using the sets of walks 
between vertices. 

The algorithm used for testing of three big groups of SRG-s, i.e. the 
groups with parameters (35,16,6,8), (36,15,6,6) and (37,18,8,9) can be 
informally described as follows: 

1. generate the sets of local complements Sq = {Gl(u),u G V(G)} for all 
graphs G in the group of nonisomorphic graphs with the same parame- 
ters 

2. compute the determinants of adjacency matrices of graphs in Sq, or- 
der the sequence of determinants by magnitude; next we partition the 
group of nonisomorphic SRG-s into the equivalence classes of graphs 
with the same (ordered) sequence of determinants of local complements; 
the graphs belonging to a class containing more than one member are 
indistinguishable by this graph invariant 

3. now it is sufficient to apply the basic variant of algorithm only for test- 
ing isomorphism of pairs of graphs, contained inside the same equiva- 
lence class. 

The group of 3854 graphs with parameters (35, 16, 6, 8) has 44 equivalence 
classes with more than one member, 42 of them have exactly two members 
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and there are two classes with four members. In the group of 32548 graphs 
with parameters (36, 15, 6, 6) we found 160 equivalence classes (152 with two 
members, 6 with three members and 2 with four members). For the last 
group of 6760 graphs with parameters (37, 18, 8, 9) there are 3379 equivalence 
classes, each with two members. The basic algorithm applied to such small 
groups of graphs is really computationally inexpensive. 



4 Conclusions 



The presented algorithm for isomorphism testing can efficiently recognize all 
known small nonisomorphis SRG-s, thus there is a hope it will perform well 
also on bigger strongly regular graphs. This needs further investigation. 
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